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Abstract 

We discuss structural and thermodynamical properties of Baxter's adhesive hard sphere model 
within a class of closures which includes the Percus-Yevick (PY) one. The common feature of all 
these closures is to have a direct correlation function vanishing beyond a certain range, each closure 
being identified by a different approximation within the original square-well region. This allows 
a common analytical solution of the Ornstein-Zernike integral equation, with the cavity function 
playing a privileged role. A careful analytical treatment of the equation of state is reported. 
Numerical comparison with Monte Carlo simulations shows that the PY approximation lies between 
simpler closures, which may yield less accurate predictions but are easily extensible to multi- 
component fluids, and more sophisticate closures which give more precise predictions but can 
hardly be extended to mixtures. In regimes typical for colloidal and protein solutions, however, 
it is found that the perturbative closures, even when limited to first-order, produce satisfactory 
results. 
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I. INTRODUCTION 



In spite of their highly ideahzed character, 'sticky hard sphere' (SHS) models continue 
to attract considerable attention in the recent literature on mesoscopic fluids. The rea- 
sons are the following. On the one hand, these models are simple enough to allow analytic 
solution of the Ornstein-Zernike (OZ) integral equation of the liquid state theory within 
some appropriate approximation. On the other hand, their potentials, albeit crude, exhibit 
realistic physical features (repulsion plus surface attraction) typical of the interactions be- 
tween particles of many macromolecular fluids. Furthermore, the SHS models have already 
been successfully employed in studies on colloidal suspensions, micelles, protein solutions, 
microemulsions, flocculation, and crystallization [liSSyi- 

The first and most known model (SHSl) was proposed by Baxter . Its pair po- 

tential adds to a hard sphere (HS) repulsion an attractive square- well (SW), which becomes 
infinitely deep and narrow according to a limiting procedure (Baxter's 'sticky limit') that 
keeps the second virial coefficient finite j^. For pure fiuids the starting potential is 
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where r is the distance between particles, cr the HS diameter, £swi and R — a denote the 
depth and width of the well, respectively, and (3 = {kBT)~^ ( being Boltzmann's constant, 
and T the temperature). The dimensionless parameter r > is an arbitrary monotonically 
increasing function of T and a is a measure of 'stickiness' between particles. In the 
limit T — s> 00, 0swi(^) becomes the HS potential with diameter R. The logarithmic form of 
^swi serves to generate a simple expression for the Boltzmann factor, e(r) = exp [— /50(r)], 



I.e. 



eswi(r) = exp [-/3<^swi('")] 



0, < r < cr 
aR {R-ay^ , a <r <R 

1, r> R 
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The 'sticky limit' amounts to the hmit R — a ^ 0, and the second virial coefficient, 
B2 = 2n dr r2 [1 - e(r)] , becomes Bf^^^ = hm^_.,+ Ef^i = (2n/3) [l - (4r)"^] . 

The OZ equation for the SHSl model admits analytical solution within the Percus-Yevick 
(PY) approximation 0, 0, However, the difficulty in applying the multi-component 

SHSl-PY solution to polydisperse colloidal mixtures with a large number p of components 
(due to the presence of p{p + l)/2 coupled quadratic equations |9| to be solved numerically) 
has led to two more recent attempts of finding an alternative SHS mode 



analyticall y tr actable even in the general multi-component case 
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and SHS3 models 



24l2a,|2g. We shall refer to them as SHS2 
j2lL |22| . In SHS2 the adhesive part of the potential is generated from an attractive Yukawa 
tail, whose amplitude and inverse range become infinite in the sticky limit, while their 
ratio remains constant. The SHS3 model has been proposed to replace SHS2, overcoming a 
recently discovered pathology while retaining some useful features of the previous model. 

As remarked in Ref. |ll|, the crucial difference between SHSl and SHS2 lies in the 
attractive part of the interaction, since the area of the SHSl well vanishes as {R—a) ln(i?— o") 
when the width R — a goes to zero, while the corresponding integral remains finite in the 
SHS2 case. 

For SHS2, the OZ equation was analytically solved within the mean spherical approxima- 



tion (MSA) 12,ll3|, and the multi-component SHS2-MSA solution is indeed much simpler 



than its SHSl-PY counterpart. 

On the other hand, Baxter's SHSl model was never solved analytically within the MSA 
or other closures different from the PY one. Note that in the MSA the main input for 
the direct correlation function is the potential SWl itself, which has, in the sticky limit, a 
weaker (logarithmic) divergence than the Boltzmann factor. Thus it would be interesting 
to investigate other approximations which consider, more correctly, the Boltzmann factor 
rather than the potential, in a way similar to the PY case. 

Motivated by this scenario, in this paper we would like to revisit Baxter's model within 
a new framework. We consider a general class of closures which, as the PY, hinge on the 
vanishing of the direct correlation function beyond a finite range R. Each component of the 
class is identified by a different approximation to the cavity function in the intermediate 
region [a, R). We shall show that this difference survives even in the limit R a and hence 
it gives rise to different solutions, each of them being characterized by a different expression 
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for a single parameter q^. 

In this class we first consider a subset of closures which is based on density expansion 
of the cavity function. All these closures yield a simpler solution than the PY one, and yet 
they provide a satisfactory description on regimes frequently encountered in experiments 
on colloids and globular proteins. The solutions corresponding to these closures could be 
straightforwardly extended to multi-component fluids, unlike the PY case. 

In addition, we identify a second subclass of closures as a possible route to improve over 
the PY approximation. 

Finally, we provide a careful discussion on the equation of state derived from this class 
of solutions. In this process, we critically discuss the correct way of taking the 'sticky limit' 
and its influence on the predicted structural and termodynamic properties of the model. A 
short-cut, but not fully equivalent, procedure may in fact generate an incorrect result for the 
virial equation of state (EOS). We point out that pitfalls of this sort have already occurred 
in the case of SHSl mixtures Such a problem should be clearly borne in mind when 

investigating thermodynamic properties of more complex fluid models which include sticky 
interaction terms in the pair potential. 

II. BASIC THEORY 
A. Integral equations 

For a one-component fluid of molecules in a volume V, with spherically symmetric 
interactions, the OZ equation is the convolution relationship 



which defines the 'direct correlation function' (DCF) c(r) in terms of the 'total correlation 
function' h{r), deviation of the 'radial distribution function' (RDF) g{r) from its ideal gas 
value of 1, i.e. h{r) = g{r) — 1. Here, the asterisk represents three-dimensional convolution. 



p = N/V is the number density, and the convolution term 7 (r) is the 'indirect correlation 
function', continuous with its first three derivatives everywhere, even when the potential has 




(4) 




Moreover, r refers to the distance between particles. 



4 



a hard-core part. The OZ equation can be cast into the Baxter form p 

j'OO 

rc (r) = —q\r) + 27rp / du q [u] g' (r + m) , (5) 

Jo 

POO 

rh{r) = —q'{r) + 2np / du q (u) {r — u) h {\r — u\) , (6) 

^0 

where the 'factor correlation function' g(r) is an auxihary function, and the prime denotes 
differentiation with respect to r. Solving the Baxter equations is tantamount to determining 
g(r), from which c (r) and h{r) can be easily calculated. To solve the OZ equation, some 
approximate 'closure' relating c(r), h{r) and the potential (r) must be added. 

B. Closures 

For hard-core fluids, whose (p may be written as the sum of a HS contribution and a 
tail, i.e., (r) = 0hs (^) + 0taii (r) (with 0hs (r) = +oo for r < a, and = for r > cr), an 
approximation is required only for r > a. In fact, in this case both c(r) and g{r) are exactly 
defined inside the core, i.e.. 



g{r) = 0, < r < 0", 

c(r) = - [1 + 7 (r)] , < r < a. 

Let us consider some known approximations to the exact DCF: 
i) the mean spherical approximation (MSA), 



(7) 
(8) 



cmsa (r) = -(3(1) (r) , r > a. 
ii) A modified mean spherical approximation (mMSA) 22 , 



27 



(9) 



CmMSA {r) = f{r), r>a. 



where / (r) = e(r) — 1 is the Mayer function, 
iii) The Percus-Yevick (PY) closure. 



(10) 



CPY(r) = /(r)[l + 7(r)], r > 0. 



(11) 
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It is well known that the MSA assumption is rather accurate at large r values, but not 
in the region just outside the hard core. For short-ranged potentials, one thus expects that 
the MSA performance could be improved by using the mMSA which replaces —(3(j) (r) with 
/ (r). Moreover, the mMSA yields the correct zero-order term in the density expansion of 
the DCF, since Cexact ij') f (r) when p ^ 0. 

Note that, when applied to Baxter's SWl potential, all these closures imply that c(r) 
vanishes beyond the range of the potential. In the present paper we will thus focus, for this 
model, on the class of all closures which can be written as 



CGPY(r) 



- [1 + 7 (r)] , < r < a, 

Cshrink (r) , a <r < R, (12) 

0, r>R. 



We shall refer to this expression as generalized PY (GPY) approximation. We remark that 
in Eq. (jl2|) the novelty lies on the fact that it does not univocally define a single closure, 
since different approximations can be obtained by different definitions of Cshrink {f)-, the DCF 
in the shrinking interval {(J^R). As it turns out, these differences persist even in the sticky 
limit R ^ a\ 

Strictly speaking, all closures of the GPY class, when different from the PY one, are mixed 
closures, being a combination of the PY approximation for r > i? with a different approxima- 
tion to c(r) for 0" < r < i?. For instance, Cshrink ij) = — /?0taii ij") yields a MSA/PY closure, 
while one recovers the PY approximation (PY/PY) when Cghrink = [1+7 We re- 
call that combinations of closures, such as hypernetted chain - mean spherical (HNC/MSA) 
and Percus-Yevick - mean spherical (PY/MSA) approximations have already been employed 



in studies on electrolytes with sticky interactions |29l . 



C. General analytic solution 

Let us apply the GPY approximation to the model defined by Baxter's SWl potential. 
Different approximations to Cshrink (t) generate, of course, different solutions (which we are 
not able to write explicitly in analytic terms). It is however surprising that, as shown in 
Appendix A, the differences among all these solutions drastically reduce in the sticky limit 
all being included in one common parameter, while the r dependence becomes identical. 
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As a consequence of this degeneracy, the OZ equation for the SHSl model admits, for all 
members of the GPY class of (mixed) closures, the following general analytic solution 



'?shsi-gpy('" ) 



-air"^ - (T^) + ha{r - a) + q^o^ 



9 (a 



for r > 0, 



(13) 



1 + 2r] 12q^T] 



31] 
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— ysHsi-GPY(o-) = lim ?/swi-gpy(o-^) 



(14) 



(15) 



12 12r — V / R^a' 

where ysHsi-GPY(c") denotes the cavity function at contact and, for simplicity, the subscript 
SHSl-GPY has been omitted in a, b, and q^j. 

n 

The GPY solution includes and extends the SHSl-PY expression found by Baxter .6], and 
formally resembles the MSA one for the SHS2 model 3]- In fact, as regards the dependence 
on r, the functional form is identical in both the mentioned cases. The distinguishing feature 
of each element of the GPY class lies, however, in the parameter dependence of qa{fli'^)i 
that mantains a 'memory' of the situation before the sticky limit is taken, depending upon 
the closure chosen for Cshrink ij) (and upon the model) (31 1. 

Note that the sticky limits of yswi-GPY(c'''') and yswi-GPY(-R''') are in general different, 
since different approximations to y{r) may be used for a < r < R and r > R, respectively. 

The function gsHsi-GPY('") is discontinuous at r = a, since gsHsi-GPY(o"^) 7^ 
?shsi-gpy(o""*') = 0, and this jump leads to 5'shsi(^) ~ ['^'" + ^^1 d {a — r) — q^a'^S^r — a). 
Different approximations correspond to different expressions for q„ = gsHSi-GPY(o" )/o"^- 

It is also of some interest to write down an explicit expression for the DCF inside the 
core, which can be derived by using Eq. (0), and reads 



cshsi-gpy('^ ) 



cr 



-12r] q; + Urj 

r 

qa cr6{r - a), 



(a + bf 



< r < a. 



aqa 



r r] 2 

a 

a 2 \a 



r\3 



9{a-r) 



(16) 



Each solution gsHsi(^) of the GPY class is fully characterized by its approximation to 
|/SHSi(o")- In the following we thus discuss some particular solutions focusing on the cavity 
function. 
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D. Cavity function 



The cavity function y{r) is defined by 

g(r) = e{r)y{r). (17) 

Sometimes it is also expressed as 

y{r) = l + -f{r) + d{r), (18) 
which defines the 'tail function' d{r) js^ls^, involved also in the exact relationship 

c(r) = /(r) [1 + 7(r)] +e(r)rf(r). (19) 

Note that the PY approximation amounts to set d{r) = Vr > 0. 

The exact y{r) is a well-defined continuous function of r with at least 2 continuous 
derivatives, even when the potential (p (^) itself has a hard-core region or other discontinuities 
[s^ . Since in Baxter's model the cavity function at the contact ?/(cr) plays a crucial role in 
characterizing a particular approximation within the GPY class, it is important to know, 
for SHSl, the first three coefficients of the expansion 

l/exact(0-) =yO + VlV + 1/2'?^ H (20) 

( rj = pvo, vq = (7r/6)cr'^ being the packing fraction and the particle volume, respectively), 
which can be derived from the SHSl exact virial coefficients Q| (see Section V.B and Eq. 
dSni)). One finds 



yo = 1, 

(21) 



y2 = 4. 591 - 4. 076 t + 1. 518 - 0.232 
+1. 056 X 10~2 t4 ^ 1. 302 5 X 10'^ 
where 

t = T-\ (22) 

Note that, within the GPY class, mixed closures imply that ?/gpy('") = 1 + 7 (r) = ypyir) 
when r > a, whereas ?/gpy(^) for a < r < R may largely differ from the PY expression. 
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A second crucial remark is that, while yexa.cti'f') is continuous everywhere, a cavity function 
derived from an approximate closure may be discontinuous. This drawback occurs, for 
instance, in the MSA (see Eq. (|27|) below) when the potential has some discontinuities 
outside the core. In the case of SW fluid, yusAir) may thus be considered continuous at 
r = a, but it is discontinuous at r = i? (a fact which has remarkable consequences on the 
derivation of the corresponding MSA virial equation of state p&l). 



E. Correct way of taking the sticky limit 

Before ending this Section, we would like to stress also a general methodological point 
related to the way the sticky limit is performed, in view of its importance in the evaluation 
of the virial EOS. The correct procedure - hereafter referred to as sticky limit at the end - 
involves three steps: first, solve the OZ equation, within some approximate closure, for the 
'starting' potential (ps'SNiir). Second, compute structural and thermodynamic properties of 
the SWl fluid. Finally, perform the sticky limit. 

In Baxter's case, an immediate difficulty arises from the lack of analytic solution for SW 
fluids. Remarkably, one can find the correct SHSl-PY solution even in the absence of a 
complete solution of the SWl-PY model, by using perturbative expansions in powers of 
R — a. The starting point is 



where the expression Q for es\Ni{r) has been inserted into Eq. (fTTj) . Note that no approx- 
imation is contained in Eq. (PHj) . if yswi(^) is exact. The use of y-pY{.f) = 1 + yields 
fi'swi-PY('")- On inserting this expression into Eq. (jHl), one can determine the necessary 
properties of the unknown gswi-PY('") (see Appendix A) and hence, in the sticky limit, of 
gsHSi-PY(r). 

However, the evaluation of structural and thermodynamic properties involves operations 
such as derivatives and integrals, which do not necessarly commute with the sticky limit. 
An underestimation of this fact has led to the existence of some incorrect results in the 
literature when the sticky limit was carried out already at the level of the Boltzmann factor 




5'swi('^) = 6{r ~ R)+ a 



R 



9{r-a) e{R-r) i/swi(r). 



(23) 



R-a 



9 



{sticky limit from the outset), i.e. 

esHsi(^) = lira esy^^(r) = e (r - a) + a a5+ (r - a) , (24) 
which imphes that [36[ 

fi'SHSi('") = 0{r -a) ysRsiir) + a aysHsilo") h (r - a) , (25) 

csRSiir) = a aysmiify) h (r - a) , for r>a (26) 

(note that in Eq. one cannot write esnsil'^) = exp [— /^^shsiI'^)], since 0shsi(^) = 

hm/j^CT 0swi(^) diverges logarithmically). This short-cut procedure of 'sticky limit from the 
outset' is rather common in the literature on systems with sticky terms in the pair potential 
(not only SHS fluids). We stress that this approach, albeit computationally simpler, is 
not always equivalent to the 'sticky limit at the end' (which remains the unique correct 
procedure) and may be dangerous in calculations of thermodynamic properties. 

In view of this scenario, we have performed an accurate re-analysis of Baxter's work and 
found that the way of taking the sticky limit does not affect q{r), as well as any quantity 
computed from it without r-derivatives, i.e., for instance, the structure factor S{k), the 
compressibility and energy EOS's, the internal energy and the Helmholtz free energy. This 
is not the case of the virial EOS as it will be shown later. 

Being aware of the danger hidden in the procedure of 'sticky limit from the outset', Baxter 
used it only whenever possible In the present paper, we follow only the surely correct 
procedure of 'sticky hmit at the end'. 

III. PARTICULAR SOLUTIONS 

In addition to the PY approximation, we now consider other GPY solutions relevant to 
mixed closures ■ ■ ■ /PY, where the dots refer to the approximation chosen for Cshrink('"), the 
DCF in the shrinking interval (o",i?). 
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A. MSA/PY or MSA closure 



Let us start from Cghrinkl'") = — /90taii For any square- well potential, and in particular 
for Baxter's SWl, this MSA/PY mixed closure coincides with the pure MSA one. In terms 
of cavity function, one can write 

Vmsa (r) = exp [/50taii (r)] [1+7 (r) - Z^^taii (r)] for r > a, (27) 

and thus 

[1 + 7SW1-MSA (r) + In f^] , a < r < R, 
1 + 7SW1-MSA (r), r> R, 

which is clearly discontinuous at r = i?. In the R ^ a limit one gets: 

2/SHSl-MSA(cr) = 0, (28) 

since the factor 1 + 7shsi-msa (o") remains finite, while the remaining term in ?/swi~msa (o"), 
involving the product of the well- width R — a times the potential at contact, vanishes due 
to the weaker (logarithmic) divergence of — /50swi (cr) (on the contrary, R — a times the 
Boltzmann factor would yield a non-zero finite result). As a consequence, one finds 

^MSA^Q^ ^ gsHsi-MSA(r) = gHS-PY(r). (29) 

This means that the SHSl-MSA factor correlation function differs strongly from its SHS2- 
MSA counterpart and coincides with the HS-PY one. In other words, for Baxter's potential, 
the MSA is a far too poor approximation, unable to include into q (r) any effect due to the 
surface adhesion. 

On the other hand, we can anticipate that the SHSl-MSA virial pressure does not coincide 
with its HS-PY counterpart, but diverges! As discussed in detail later on (Section V.C and 
Appendix B), this is due to the discontinuity at r = i? in the cavity function yswi-MSA('") 
within the MSA approximation. Hence a particularly careful treatment of the jumpin 
the cavity function, appearing in the virial equation of state, is required in this case j35|, 
unlike all those closures which preserve the continuity of y{r) (e.g. PY). As a result, we 
find the somewhat paradoxical situation of having the limit of the factor correlation function 
?shsi-msa('") identical with the hard sphere one qHS-Pvi^), whereas the corresponding limits 
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yswi-MSA [r) 



of the virial equations of state are different. This unexpected result may appear puzzhng at 
first sight. However we attribuite it to the strong thermodynamic inconsistency generated 
by the MSA in the sticky hmit. 



B. mMSA closure or CO-approximation 

Henceforth the mMSA will be denoted as CO-approximation, to emphasize that it cor- 
responds to the exact zero-order approximation, Cghrinkl'") = f {f')^ for the DCF in the 
shrinking interval outside the core. 

Now yco (r) = 1 + 7 (r) exp [/?0taii (r)] for r > a, and 

1/swi-co (r) = 1 + ^ 7SW1-C0 (r) , a < r < R. 
an 

This implies that 

2/SHSi-co(o") = 1, (30) 

e = ^ , (31) 

which depends on temperature, but is density-independent. 

It is interesting to remark that the resulting 5'shsi-co(^) solution coincides with the 



mMSA solution of the SHS3 model 



22j . provided that the reduced temperature T* of that 



potential is replaced with Baxter's r. It is also wort 



1 noting that the SHS3-mMSA solution 
is formally coincident with the SHS2-MSA one 

[211, Q- As a consequence of this mapping, 
most of the MSA results obtained in the past for SHS2 (which turns out to be ill-defined 
rom a thermodynamical point of view) are found to still hold within a correct framework 



221. 



C. Cl-approximation 

The most natural extension of the previous approximation is the inclusion of the next 
order correction in the density expansion. We shall denote this as Cl-approximation, which 
is defined by 

Cshrink(r) =/(r) [1 + 71 (r)p], (32) 
12 



where 71 (r) = (/ * /) (r) is the first-order coefficient in the density expansion of 7 {r,p). 
This approximation is exact up to first-order in density. As a consequence, we find 

l/SHSi-ci(cr) = 1 + yiit) T], (33) 
with Ui being the exact first-order coefficient given by Eq. (j2U). The corresponding q^^ 
follows immediately. 

Clearly, higher orders terms in the density expansion of Cshrinkl'^) could be considered. 
However, even the inclusion of the next order correction would require a much more elaborate 
analysis compared to the first-order one. Furthermore it is reasonable to expect that higher 
order, t-dependent, terms should be introduced in a way balanced with the HS part, in view 
of the slow convergence of the perturbative series. For these reasons, in this paper, our 
analysis will be limited to the ffist correction only. 



D. Baxter's PY solution 



Next we consider the PY closure, which cannot be deduced from CO by adding a finite 
number of perturbative terms. 

In this case the approximation to y{r) is 'continuous', in the sense that the same closure 
is employed both in (a, R) and {R, +00). Eq. (fT^ then becomes 



t 



(34) 



where the contact value of the PY cavity function can be determined analytically by solving 
a quadratic equation, usually written in terms of A. Only the smaller of the two real solutions 
(when they exist) is physically significant Q]. The result reads 



?/SHSi-PYlcrj 



1 + T^^ + 



(35) 



(36) 



where yo^cr) is the HS-PY result 

l/HS-PYv(^) = (l+W2)/(l-r/)^ 

It is interesting to compare the ffist terms of the expansion of ysHsi-PY(<7) in powers of t] with 
the exact results given by Eq. PT|) . One finds that y^^^^'^^ = 1 
while 



„ exact SHSl— PY exact 
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SHSl-PY ^ 4 _ !^ + 11^2 _ 1 ^3 1^4 

^2 2 12 4 72 

= 4.0 - 3. 5 t + 1. 417 - 0.25 t^' + 1. 389 x 10"^ t^. (37) 

Note that in y^^'^ there is a t^-term, missing here. 
E. Beyond the PY closure 

By exploiting the Baxter-OZ equations, one can evaluate 1 + 7shsi-gpy(c") and find (see 
Appendix B) 

7SHS1-GPy(0") = ?/HS-PYv(cr) " 1 " Y~^^ ^ 12^^^' ^^^^ 

An improvement over the PY approximation can be obtained by recalling Eq. (fTH|) and 
writing 

1/shsi-gpy(o") = 1 + 7shsi--gpy(o") + c?(cr) 

V 1 2 

= yHS-PYv(cr) - ^ ySHSl-GPY(cr) + [t ysHSl-GPY(o")] + rf(o-)(39) 

We stress that in Eq. ()39|) the presence of d{(j) is our key for improving over the PY closure. 
If we are able to give an explicit approximation for d{cr), then ?/shsi-gpy(c'") can be computed 
(analytically or numerically) from Eq. ()39|). 

In the following we shall attempt to motivate a particular form of d{a) using a combination 
of exact heuristic arguments. Our guidance in this task hinges on two observations. First, 
we note that density expansion of Eq. (fTHj) suggests that d^cr) starts with a p^-term, and 
hence it is negligible at low density. The PY case indeed corresponds to set d{a) = 0. 

On the other hand, it is well known that the PY approximation yields a description of the 
cavity function for HS which is far from being fully accurate. This can be quickly established 
from the density expansion of Eq. which yields y^^~^^^ = 4, markedly different from 
^HS-exact _ 4 59]^ (^ggg (pT|l ). In order to cope with both the above remarks, we propose 
the following simple expression for d{a), which will be referred to as GPYl, 



14 



Substitution in Eq. yields Eq. again, but now with 

2/o(^) = ycsicr) - {t/2 - e/l2) r,\ (41) 
ycs{o) = (1-W2) / {l-ri)\ (42) 

The HS part of (iGPYi(c") has been chosen in order to get jj/nsftr), the well-known Carnahan- 
Starling expression, which is more accurate than ?/hs-pYv (f^) |22l ■ On the other hand, the 
i-dependent part of (iGPYi(c) ensures that now the r;— expansion of ysHsi-GPYi(c") yields 

^SHSl-GPYl ^ 4 5 _ 4 ^ ^ 15 ^2 _ 1^3 ^ 1 ^4^ (43) 

where the first three terms are closer to the exact ones than their PY counterparts. 

In spite of the fact that the approximation given in Eq. (jin|) is not fully physically 
motivated, it nevertheless does represent an improvement over the PY one, as it will be 
shown in Section VI where numerical tests are reported. 

IV. STRUCTURE 

The structure of the fluid is described by g{r), or by the structure factor, 

S{k) = 1 + ph{k) = [l- pc{k)]-^ = [g {-k) Q (A;)] , (44) 

where k is the magnitude of the scattering wave vector, h{k) and c{k) denote the three- 
dimensional Fourier transform of h (r) and c(r), respectively, while Q{k) = 1— 27rp q{k) , with 
g (fc) being the unidimensional Fourier transform of q (r). Using Eq. (fT^ . the calculation of 
S{k) is straightforward and yields 



where 



SsRsi-GPY (ka) = ^ m2^mo - fu 

[1 - 12r] Qcos {kcr)\ + [127] gsin {ka)\ 



gcos (x) = h (x) +h h (x) + go h {x) , 

gsin {x) = J2 (x) + 6 Ji (x) + go -^0 {x) , (46) 



go = g(0)/a2 = - [2(l-r/)]-^ + g^, 
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(47) 



and 



du u'^ COs{xu) , Jn{x) 



du u" sm{xu). 



(48) 



V. THERMODYNAMICS 



The thermodynamic properties can be obtained through the usual compressibihty, energy 
and virial routes to the EOS. In any approximate theory, however, these results do not agree 
with each other, due to thermodynamic inconsistency. 

The exact values of the first virial coefficients in the //-expansion of the SHSl pressure 



are 



/3Pvo = 7] + B2rf + B:,rf^ + B^r]'' + 



(49) 



is'' ' 



(50) 



S2 = 4 - t, 
^3 = 10 - 5t + - 

Bi = 18.36477 - 13.77358 t + 6.114 - 1.518 

+0.17398 - 6.33514 x 10"^ - 6.51271 x 10"^ 
In the following we will discuss results for the MSA, CO and CI closures, and compare 
them with the PY ones. Thermodynamics of the GPYl approximation will not be considered 
in this paper. 



A. Compressibility route 



In general, the compressibility (C) equation reads 

where the factorization given in Eq. ()44|) as well as the definition of a have been exploited. 
Integration with respect to density then yields the compressibility EOS, i.e. for Z = (3P/p 



(51) 



drj. 



(52) 
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For our model and within the GPY class of closures, Eqs. ()14|) yields 



asHsi-GPY = — TT 7, Kv, t) + ) A (r/, t). (53) 

(l-r/) (1-??) 

In particular, for the MSA one has the simple result 



7SHSI-MSA 7HS-PY 1+77 + 77^ , s 

=Zc = (i_^)3 • (54) 



which is consistent with that reported in Eq. 

Then, for the cq- approximation (mMSA) compressibility EOS one finds 

^SHSi-co ^ ^ ;5pSHSi-co y^ = t^r]+ (4t + 2t') ln(l - r/) 



+ (1 + 4t + t') ^ 



-3(1-0(^)^3(^1. (55) 



while the ci- approximation gives a more complicated expression 



7SHS1-C1 ^ _ 7SHS1-C0 ^ 
V - V 

+ (4 + 4t + 3Ai) \iri + {t + Xi) W + -Xi^V^ 
+ (14 + 6t + 4Ai) Ai In (1-77) 

+ (10 + 2t + Ai) Ai-^ - 3Ai (j^) , (56) 
where Ai = t yi{t). These results can be compared with the well known PY ones 

^SHSl-PY ^ ^HS-PY _ ^ ^ + ^3 _ (57) 

Explicit numerical tests will be reported in Section VI. 

It is worthwhile noticing that use of Eq. (j3(J|) into Eq. (|57|) does not yield back Eq. (|55|1. 



B. Energy route 

The Helmholtz free energy A can be computed from the cavity function at contact. In 
fact, following Baxter [3], it can be shown that 

/3 (AsHsi - ^Hs) /N=-r] [ i/sHsi((x) dt. (58) 

Jo 
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We stress that this relationship is exact, independent of any closure. Once that A is known, 
the energy EOS can be determined via Ze = f] d{(3A/N)/ Or]. 

Inserting the analogue of the ?7-expansion Eq. (pUj) for ysHSi(c") into Eq. (fHH|). one gets 
the exact relationship between virial coefficients and yn{t), namely 



?HS 



Bn -B^^ = -{n-l) du yn-2{u), n > 2, (59) 

which we have exploited to calculate the exact first yn{t)^s from the exact first -B„'s given 
by Post and Glandt S^- 

/^From the contact values calculated in Section III.B one finds 

/N = PAns/N (60) 

/5^SHsi-co/iV = PA^s/N - t T], (61) 

/?AsHsi-ci/iV = PAns/N -trj+ (^-5t + - j^t'^ ^ (62) 
and, the corresponding expression for Ze are: 

^SHS3-C0 = Zhs - t V (63) 



^SHS3-ci = z^s-tv+ [-5t + e- J (64) 

Within the PY approximation, analytic expressions for Ashs3-py and Z^ ~ are also 
available [l^ . 

It is worth noting that the CO-energy results are mean-field ones: Z^ ~ predicts only 
the second virial coefficient correctly, whereas its higher-order virial coefficients coincide 
with the HS ones. On the other hand, the CI and PY closures give also the exact i?3. 

C. Virial route 

The determination of the virial (\^) EOS 

27r /"°° 

Zv = l + —pj^ dr r'g{r) [-/30'(r)] . (65) 
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is the most delicate one. The PY result is 



7SHS1-PY 



V 



l + 2r/ + 3r/2-^(l + 5r7)Ao + ^A2 



(66) 



3r(l -r/)3 

with the shorthand Aq = Xpyfj^l — f]). 

As remarked in Section II. E, a special care must be exerted in evaluating the sticky limit 
of this case. As no explicit derivation of Eq. ()fi6j) has been reported in the literature, to the 
best of our knowledge, we shall discuss this point in some detail here. 

The starting point is the exact expression , 



Zf SI = lim zswi = l + 4r]\ ysnsiicr) - a [3ysHsi(^) + ^ Hm y'swii^)] } ■ (67) 

R^a K. L R-^a J J 

The evaluation of the last term is the crucial point. To our knowledge. Bravo Yuste and 



Santos 



41| were the first to clearly stress that, in general, ImiR^a UswA^) cannot be 



replaced with ?/shsi('^) = [lini-R^o- l/swi(^)]r-=o- ; analogous quantity corresponding to the 
procedure of 'sticky limit from the outset', since sticky limit and r— differentiation do not 
commute (see Appendix B). In fact, on replacing the last term in Eq. (jHIj) with 1/shsi(^)' 
one obtains an incorrect virial EOS, which differs from Eq. ()66|) (in particular, the Ajj-term 
is missing). 

In some sense, a 'memory' of the original potential is retained even after the sticky limit 
has been performed. The impossibility of getting the correct derivative from y{r) once the 
sticky limit has been taken was also recognized by Seaton and Glandt 

In short, the correct PY virial EOS cannot he obtained if the sticky limit is taken 'from 
the outset'. On the other hand, no such drawbacks occur when evaluating other structural 
or thermodynamic properties in the absence of r— differentiation. This important point 
has sometimes been overlooked in the past literature. For instance, we mention here the 
extension of Eq. to mixtures as proposed by Barboy and Tenne Q|. Their multi- 

component PY virial equation of state, based upon the 'sticky limit from the outset', is not 
correct, since it does not reduce to Eq. ()66p in the one-component limit, as it should. A 
similar pitfall appears in the virial pressure ofRef. Q. 
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In Appendix B it is also shown that 



^SHSl-MSA ^ ^SWl-MSA ^ (gg) 

R^cr 

i.e., in the sticky hmit the MSA virial EOS differs from Zy^~^^ and diverges. 

Although this result may appear odd at first sight, it is nevertheless consistent with our 
previous discussion. The virial EOS depends not only on g(r) after the sticky limit, but also 
on the 'tail' of the original potential. 

It is also worth stressing that the evaluation of the MSA virial EOS for SW fluids embodies 



a further sub 

2/swi-msa('^) 



;le point not present in the PY virial EOS, because of the discontinuity of 
(see Appendix B). 

One could reasonably expect the above divergence to be a consequence of the poorness 
of the MSA closure, which should disappear by resorting to more accurate approximations. 
Unfortunately, the virial EOS diverges within the CO-approximation too (see Appendix B): 



^SHSl-CO ^ ^SWl-CO ^ (gg) 

_R— XT 



and thus higher-order closures should be considered. 



VI. NUMERICAL RESULTS 



In order to give a flavor of the difference appearing in the different approximations, we 
present here some results and an interesting comparison with Monte Carlo (MC) simulation 



data 
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A. Closures simpler than PY 

Figures 1 and 2 depict MSA, CO, CI and PY results for the structure factor S{k), com- 
pared with the corresponding MC data |3 . 4^ . Four thermodynamic states are considered 
and we have displayed the corresponding results in order of decreasing r. The simulation 
data for S{k) stem from a recent improved re-analysis 143] of the calculations by Kranendonk 
and Frenkel These data are less susceptible to finite size effects than the original ones 



42|. 
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As MSA curves coincide with the HS ones, one can clearly appreciate the effects of 
.Ace a..ac.io„ p.e.n. ,„ .he SHSl .ode,. A. . .now„ fl [j at M,. valuer of . 

repulsive forces predominate and the structure factor of sticky hard spheres resembles that 
of hard spheres, although heights and positions of peaks are quite different (see Figure 
1). Lowering r (i.e. decreasing temperature or increasing stickiness) leads to a different 
behaviour characterized by a strong increase of S{k = 0), which is related to the isothermal 
compressibility and the density fluctuations. The observed development of such a new peak 
at k = (see Figure 2) signals the approach to a gas-liquid phase transition (S'(O) diverges 
at the spinodal curve). While there is a good overall agreement between MC data and PY 
predictions, one can observe that the CO and CI curves lie rather close to the PY ones. In 
addition, in Figures la), lb) and 2c) one clearly sees an improvement in going from the CO to 
the CI approximation, whereas in Fig. 2d) the CO results are nearly coincident with the PY 
ones, at least in the considered interval of parameters. In Figure 2 the largest discrepancies 
among the different approximations occur near the origin ( see 2c) ). 

Figures 3 and 4 illustrate the behaviour of the compressibility pressure, /SPccr^- In all the 
considered cases, the performance of the CO and CI approximations is quite good (again, CI 
is slightly better than CO). The CO, CI and PY results are nearly coincident in the interval 
< rj < 0.15, which covers the concentration range typical of many colloidal or protein 
solutions. 

Figures 5 and 6 show the predictions for pressure obtained from the energy route, (3Pect^- 
Now, there are much larger differences with respect to the PY and MC results. It is 
somewhat surprising that, unlike the PY case |4^, within the CO and CI closures the 
compressibility route appears to be superior than the energy one. 



B. Closures beyond PY 

As a final illustrative example, we report here a sample of GPYl results. In Figure 7 
GPYl and PY predictions for S{k) are compared with the MC simulation for the state 
(r, r/) = (0.2,0.32). In general the PY theory tends to underestimate the amplitude of the 
oscillations in S{k). Now the improvement present in 1/2^^^ is reflected in the behaviour of 
S{k) near the origin. In this region, while the PY curve lies below the MC data, the GPYl 
one has a more correct trend. On the other hand, since the PY critical point evaluated via 
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the compressibility route ( i.e. via S{0) ) is known to be underestimated with respect to 



the exact one 



4J|, one can reasonably expect that 5'shsi-py(0) may be underestimated too. 



We have not considered thermodynamics corresponding to the GPYl closure, since such 
a task would probably require analytical/numerical calculations which go beyond the scope 
of the present paper. One could expect, however, results not significantly different from the 
PY ones. 



VII. CONCLUSIONS 

We have shown that Baxter's model (hard spheres with an infinitely deep and narrow 
attractive square-well) is analytically solvable not only within the PY approximation but 
also within a general class of closures. We have denoted this class as 'generalized PY' 
approximation, since all its elements are 'mixed closures', which coincide with the PY one 
beyond R, the upper boundary of the well, but may employ a different approximation within 
the well region (a, R). 

It has been argued that, even in the limit of vanishing width of the well, there is a 
'memory' of the approximation to y{r) chosen within the well, and this has been identified 
as the distinguishing feature of each closure of the GPY class. 

In particular, we have shown that: i) the SHSl model can be solved analytically within 
the MSA closure, but its SHSl-MSA solution has nothing to do with the SHS2-MSA one 
and is equal to gHS-PY('^); h) an interesting analytic solution can be obtained also within 
the mMSA (or CO-approximation). The resulting factor correlation function q'sHsi-co('") is 
formally identical with qsHs.s- piiu s a (r) (if r is replaced with the reduced temperature T*), 



obtained for the SHS3 model l21 
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Although the solution q{r) of the OZ integral equation, and all properties which can 
immediately derived from it, can be carried out within the unified GPY treatment, the cal- 
culation of other properties may require an analysis starting from the original SW potential, 
i.e. before the sticky limit. A detailed analysis of the non-commutativity of the sticky limit 
with other operations in some instances, notably for the virial EOS, has been provided. In 
principle, one could extend the above remarks from the SHS case to any general model with 
sticky terms in the pair potential. Within this more general viewpoint, the presence of Dirac 
(5-functions in any of the functions f{r), c(r), g{r) could be regarded as an indication that 



22 



the procedure of 'sticky limit from the outset' is imphcitly adopted. The experience of the 
present paper on SHS systems suggests that special care should be exerted in evaluating 
properties where r-differentiation is required, as in the virial EOS, since the corresponding 
results might be altered by the non-commutativity of the sticky limit. 

The GPY class of closures may be roughly split into two subclasses having the PY 
approximation as a separating borderline. The first subclass is originated by the MSA, 
which is commonly believed as one of the simplest possible closures in view of its mean- 
field character. The main drawback of this closure lies in its linear relationship with the 
potential, which prevents a correct zero-density limit of the correlation functions. We have 
cured this by modifying the closure so that the CO-direct correlation function is related to 
the Mayer function rather than to the potential, and in this case the correct zero-density 
limit is recovered. Finally, an explicit density-dependence can be perturbatively included 
in the closure, order by order, through a density expansion. We have discussed in detail 
the analytical results stemming from the zero (CO) and first (CI) order terms in density 
expansion and a comparison with Monte Carlo simulations shows that, in spite of the crudity 
of the approximations, the resulting structure factors are rather close to the MC data for a 
wide range of the parameters. 

Our main interest in the closures of this first subclass is that they are simple enough to 
allow an easy extension to mixtures and application to polydisperse fluids, unlike the PY 
case. The main disadvantage is, on the other hand, represented by strong thermodynamic 
inconsistencies, which cannot be coped with in a simple way. 

The closures of the second subclass have in fact been devised to the aim of reducing 
these inconsistencies, by exploiting the additional contribution d{a) to y{cr)- We have given 
only an illustrative example (called GPYl) to show how small differences in the contact 
value of the cavity function can produce, for instance, a different behaviour of S{k) near 
the origin. On the other hand, the exact DCF is expected to possess also a 'tail' for r > R, 
which all GPY closures neglect. The presence of this tail (which could be approximated, 
for instance, by a Yukawa function) implies that the exact solution gexact('") cannot coincide 
with gGPY(^)- These slight differences can be expected to be essential in order to achieve 
thermodynamic consistency, so that the use of the mentioned d{a) is useful but probably 
not sufficient. In this framework, it would be interesting to compare our (i(cr)-based closures 

nn 

with the completely different approach by Santos et al. p|, |41||, who proposed, for SHSl, to 
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go beyond the PY closure through a rational function approximation. 

Finally, we wish to mention that work along the lines laid down in the present paper 
is in progress to investigate, by the mult i- component extension of the CO-approximation, 
stability boundaries and percolation threshold of polydisperse colloidal fluids. 
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APPENDIX A: METHOD OF SOLUTION 

The condition c(r) = when r > R implies that q{r) = for r > i?, and q{R) = 0. Thus, 
Eq. (jni) may be cast in the form 



where the definitions G{r) = rg{r), a = 1 — 27ip du q{u), ha = 27rp du uq{u), as well 
as the core condition, Eq. (jHJ, have been employed. Here, mm{R, r — a) denotes the smaller 
between R and r — a. 

iFiom. Eqs. ()A1|) . one finds that q{r) splits into three parts: qi{r) for < r < i? — cr, 
q2{'r) for R — a < r < a, and q3{r) for a < r < R. The function q{r) is continuous at 
r = R — a, a, and R, while q2ir) is always represented by the second degree polynomial 



As -R — 0" ^ 0, the intervals (0, R — a) and (cr, R) vanish, and q{r) reduces to the sticky 
limit of q2{r). To obtain asnsi = lini_R^o- a, as well as &shsi, it is indifferent to take the sticky 
limit 'at the end' or 'from the outset' (see Section II. E), since such a limit commutes with 
the u-integration included in the definition of a and b. Because of the continuity condition 
5'2(cr~) = g3(cr+), the sticky limit of g(cr) can be evaluated from g3(r). Although the explicit 
solution is not known, integrating Eq. ()A1|) yields a formally exact expression of qsir), from 
which we get 



G (r) + g'(r) = ar + ba + 2-np 




(Al) 



\a{r^ -a^) + b<y{r -a) + q{a) [10]. 
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g3(fT+; R) = la{a^ - R^) + ba{a-R)+ R) - 2-Kp F2(a; R), 

R)= dt G'(t) = Gia+) (R - a) + -G"((T+) (i? - a) V ■ • • , 

F^[a- R) = / ^'^"^ ~ ^^^^^^ " ^?i(0+)G((T+) iR-af + ---. 

It is reasonable to assume that a and 6 tend to finite values, and the same occurs for qi{0~^), 
whereas G(cr+) and G"(a+) may diverge. Consequently, gsHSi(c") is determined by 

l-R 

lim q^{a+]R)= lim / dt G{t) = lim G^kn^],{a^]R){R- a), 

since all remaining terms vanish as i? cr (G'shrink(o"^; -R) could also be replaced with 
Cshrink(c'"'^; -R), since the difference becomes zero in the sticky limit). Using Eq. (fTTj) . along 
with Eq. Q, then yields Gswi{a^] R) = a a yswi{<^^',R) "rz^; and Eq. (fT3j) follows 
immediately. That simple general relationship allows to determine q{r) of the SHSl model 
for any closure with c(r) = when r > R. 



APPENDIX B: THE VIRIAL EOS 

1) Exact case and PY approximation. 

For any SW potential, the virial theorem is awkward since both g{r) and 0(r) are 
discontinuous at r = a and r = R. Usually, one rewrites g{r) [—(3d(l){r)/dr] as y{r) de{r)/dr 
in Eq. and uses 

^^^^ = 6{r - a) + (1 - e^^)5(r - R), 
dr 

with —e being the well depth. If y[r) is continuous at both r = a and r = R (as normally 
occurs), then one obtains 



Z^w = 1 + _p [eP^a' yswicT) + (1 - e^^R'yswiR)] (Bl) 
= l + 4r] {gsw{cT+) + {R/crf [gsw - 9sw{R-)] } , 

which holds for any SW fluid and contains no approximations. For the particular SW of 
Baxter's model, the previous exact EOS becomes 
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As the exact (as well as the PY) y{r) and y'{r) are finite and continuous everywhere, one 
can expand r^yswi{r) about r = a 



r^yswiir) = cr'^yswiicr) + [r^yswi(r)]Ua (r - a) + C((r - af). 
s expresi 

HQ 



Evaluating thi s expre ssion at r = i?, inserting it into Zy^^, and finally taking the sticky 
limit, one gets 



Zfsi = \imZ^'^' = l+'^p\a'ysusi{cr)-aa lim[rW(r)]U| 

= l + 4r]\{l-3a)ysnsi{(^)-aa lim [yswi(r)]U| • (B2) 

Once more, this result is exact, independent of any closure. 

In contrast to what happens for one-dimensional systems of sticky hard rods [4^, the 
order of sticky limit and differentiation is essential in the three-dimensional case, since these 
operators do not commute now, i.e.. 



lim ( lim ) yswiir) ^ ( lim ) lim yswiir). 

R-KT \r^a dr J \r^a dr J R^a 



This is true at least in the PY approximation, where one finds that 



lim [?/swi-PY(r)]U<x = 2/shsi-py(^) - 7^^^^^- (B3) 

To highlight the source of such a difference, we give a detailed demonstration of this result. 
Let us define Y{r) = ry{r), and consider its PY approximation, YpY(r) = r [1 + 7(r)] . Since 
7 (r) equals h{r) — c (r) , this function can be obtained by subtracting Eq. ^ from Eq. 
After expanding all integrals dt . . . in powers of R — a and using the relationship 
= (a + 6) cr — G{a'^), one takes the sticky limit and obtains 

>SHSi-PY(a) = (a + 6) a + 12r] aYia) g(0) /(T^ (B4) 



lim FsVi-pyM = « - 12r7 {a + b) q{0) - aY (a) \im q's^,_py{0) /a^ (B5) 



R- 



R- 
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(for simplicity, the subscript SHSl-PY has been omitted in both right hand sides). 

To evaluate lim/j^o- 5'swi-py(0)' '^^^ S^^^ ?swi-py(0) from Eq. ()A1|) for < r < _R — cr. 
After expanding dt q (t) G (t) up to the {R — cr)^-term (again with the help of q'{ci^) = 
{a + b) a — G{a^) ), the sticky limit can be performed and yields 

hm g^wi-PY(O) = (^6shsi-py - ^• 
On the other hand, if one takes the derivative of q'sHsi-PY(?") directly, then 

9sHS1-Py(0) = ^SHSl-PY 

In conclusion, the difference in Eq. f)B3|) stems from the following one: 



lim gswi-PY(O) = gsHsi-PY(O) - TTT^^^^- (B6) 

2) MSA closure. 

The derivation of the MSA virial EOS is not plain and thus deserves some remarks. In 
fact, one might believe that the SW-MSA virial EOS can be simply obtained by inserting 
fl'sw-MSA (R^) — fi'sw-MSA(-R^) = — /9e into Eq. ()B1|) . which is exact for a generic SW fluid. 
Unfortunately, this opinion is incorrect. In fact, as pointed out by Smith et al. ^3^], in the 
SW-MSA case the procedure used for Eq. (jBljl becomes meaningless, because ?/msa('") given 
by Eq. ipTjl is continuous at r = a, but discontinuous at r = i? where the SW potential is 
discontinuous (on the contrary, ?/exact('^) is continuous everywhere). Consequently, Eq. ()B1|) 
cannot be valid in the MSA. 

Using the relationship 



B 

dx F{x)S{x — xq) = 
A I 0, otherwise. 



i [F (xo ) + F (x+)] , if A<xo<B, 
0, otherwise, 
which extends the usual definition of Dirac's S to F-functions discontinuous at Xq, one can 



demonstrate that the proper MSA expression is 
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^SW-MSA ^ l+4,^|^s^_^,s^(a+)-/3£ {R/af ^ [(7sW-MSA(i?-) + ^?SW-MSA(i?+)] } 

= 1 + A7]\l + 7sw-msa(cx) +(3e-f3e {R/a 



1 +7sw-msa(^) + ]j^(3e 



(B7) 
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For Baxter's model, the sticky limit of tends to — oo, as — {(^Eswif — 

3) CO- approximation. 

The CO-calculations are similar to the MSA ones. Now the resulting EOS is 



7SW-C0 



= 1 + 4r7 {1 + 7sw-co(cr) + e'^' - 1 



-(3e {R/af 



1 + 1sw-cq{R) + ^e"' 



(B8) 



Unfortunately, also the sticky limit of Z^^ is infinite, with the most divergent term 
being proportional to — /Ssswic'^^^'^^ = —a ;^z^ln {a ■^^)- 
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FIG. 1: Structure factor S{ka) versus dimensionless wavector ka for various closures ranging from 
MSA to PY, and comparison with Monte Carlo results (MC) from Refs. [43,44]. a) r = 1.0 and 
rj = 0.5; b) r = 0.5 and rj = 0.4 [ here, as well as in Fig. 2, the curves of the upper part are shifted 
upwards by 2.5 units ]. 

FIG. 2: Same as in Fig. 1 with different temperatures and packing fractions, c) r = 0.2 and 
r? = 0.32; d) r = 0.1 and r) = 0.14. 

FIG. 4: Same as in Fig. 3 with different temperatures c) r = 0.2; d) r = 0.1. 



FIG. 5: Energy pressure PPeO'^ as a function of the packing fraction r] for two different values of 
temperature (a) r = 1.0; b) r = 0.5) and various closures. 



FIG. 6: Same as in Fig. 5 with different temperatures: c) r = 0.2; d) r = 0.1 . 



FIG. 7: Structure factor S{ka) versus dimensionless wavector /ccr for r = 0.2 and rj = 0.32. The 
performance of the GPYl is tested against MC data of Ref. [43] and PY results. 



FIG. 3: Compressibility pressure f3Pccr'^ as a function of the packing fraction rj for two different 
values of temperature (a) r = 1.0; b) r = 0.5) and various closures [ here, as well as in all following 
figures for pressure the curves of the upper part are shifted upwards by 2 units ]. 
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